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ABSTRACT 

Gas in galactic disks is collected by gravitational instabilities into giant atomic-molecular complexes, 
but only the inner, molecular parts of these structures are able to collapse to form stars. Determining 
what controls the ratio of atomic to molecular hydrogen in complexes is therefore a significant problem 
in star formation and galactic evolution. In this paper we use the model of H2 formation, dissociation, 
and shielding developed in the previous paper in this series to make theoretical predictions for atomic 
to molecular ratios as a function of galactic properties. We find that the molecular fraction in a 
galaxy is determined primarily by its column density and secondarily by its metallicity, and is to good 
approximation independent of the strength of the interstellar radiation field. We show that the column 
of atomic hydrogen required to shield a molecular region against dissociation is ~ 10 Mq pc~^ at 
solar metallicity. We compare our model to data from recent surveys of the Milky Way and of nearby 
galaxies, and show that the both the primary dependence of molecular fraction on column density 
and the secondary dependence on metallicity that we predict are in good agreement with observed 
galaxy properties. 

Subject headings: galaxies: ISM — ISM: clouds — ISM: molecules — ISM: structure — molecular 
processes 



1. INTRODUCTION 

The formation of molecular hydrogen is a critical step 
in the transformation of interstellar gas into new stars. 
The neutral atomic interstellar medium (ISM) in galax- 
ies is generally segregated int o cold clouds em bedded 
in a warm inter-clo ud medium (jMcKee fc Ostrikeij|l977t 
iWolfire et al1l2f)f)l . and the inner parts of some of these 
cold atomic clouds harbor regions where the gas is well- 
shielded against dissociation by the interstellar radiation 
field (ISRF). In these regions molecules form, and once 
they do star formation follows. 

A full theory of star formation requires as one of its 
components a method for expressing in terms of observ- 
ables the fractio n of a galaxy's ISM t hat is in the molec- 
ular phase (e.g. iKrumholz fc McKee 2005). No models 
published to date satisfy this requirement, but observa- 
tions have yielded a number of empirical rules for galax- 
ies' molecular co ntent. Based on Hi and CO mapping 
of n earby g alaxies 'Wong fc Bhtz (2002", hereafter WB02) 
and|Blitz_fc Rosolowskv (2004, 2006, hereafter BR04 and 
BR06) infer that the molecular to atomic surface density 
ratio i?H2 = in a galaxy varies with the inter- 

stellar pressure P needed for hydrostatic balance in the 
ISM as i?H2 cx; P^-^^^ g^j^^ thai the atomic surface density 



saturates at a maximum value of ^ 10 Mq pc~^. The ob- 
served saturation and a similar dependence of molecular 
fraction on pressure are also seen in newer surveys such as 
the HERA CO-Line Extragalactic Survey (HERACLES) 
that cover a broad er range of galaxy properties at higher 
spatial resolution (| Walter et al.l 120081 : iBigiel et all 120081 
iLerov et al.ll2008l . hereafter L08). However the physical 
origin of these patterns is unclear. The samples on which 
they are based are composed solely of nearby galaxies 
with a limited range of properties, and in the absence of 
a physical model it is uncertain how far they can safely 
be extrapolated to regimes of metallicity, surface density, 
or other properties not represented in samples of nearby 
galaxies. 

Theoretical treatments of the problem to date do not 
yet make such an extrapolation possible. A number of 
authors have considered the microphysics of H2 forma- 
tion and the structure of photodissociat ion regi ons in 
varying levels of detail (e.g. Ivan Dishoeck fc Black HM 



Black fc van Dishoeck|[l 987: Sternberg 1988: Elmegrceu 
19891: iDraine fc Bert oldi 1996: Ncufcld fc Spaans 1996 



Spaans fc Neufeldl 1997: HoUcnbach fc Ticlcns 199j; 



Liszt fc Luc as 2000; iLisztl 120021: IBrowning et al.l 12003 



Allen et al.i,2004) . but none of these treatments address 
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the pr oblem of ato mic t o mole c ular ra tios on galactic 
scales. IWvsd (|1986f ) and IWand (|1990al |bh present mod- 
els for cloud formation in galactic disks, but these both 
rely on prescriptions for the rate of conversion of atomic 
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to molecular gas that are based either on rates of cloud 
collisions or on Schmidt laws, not on phys i cal rn odels of 
H2 formation and dissociation. lElmegreenI ()1993f ) gives a 
theory of the molecular fraction in galaxies that does in- 
clude a treatment of the H2 formation and self-shielding. 
However, his model neglects dust shielding, an order 
unity effect, and it also requires knowledge of a galaxy's 
ISRF strength, which cannot easily be determined obser- 
vationally, and its interstellar pressure, which can only 
be inferred indirectly based on arguments about hydro- 
static balance. This makes the model difficult to test 
or to apply a s par t of a larger theory of star forma- 
tion. ISchavd ()2004f ) considers the conditions necessary 
to form a cold atomic phase of the ISM. The existence of 
such a phase is a necessary but not sufficient condition 
for molecule formation, so while Schaye's model provides 
a minimum condition for star formation, in makes no 
statements about what fraction of the ISM goes into the 
molecular phase able to form stars, and thus no state- 
ment on the star formation rate that is achieved once 
the minimum condition is met. 

Numerical m o dels are in a similar si t uation . 
iHidaka fc SofiI3 (|200^ and iPehipessv et"all (l2006h 
simulate galaxies using subgrid models for H 2 form ation 
similar to those presented by lElmegreenI (|1993l ). and 
show that they can reproduce so me qualitative features 
of the H2 distribution in galaxies. [Robertson fc Kravtsovl 
(|2008[ ) show that a simulation of a galaxy's ISM that 
includes radiative heating and cooling in the ionized 
and atomic phases, coupled with an approximate 
treatment of II2 formation on grains and dissociation 
by the ISRF, can reproduce the observed molecular 
content of galaxies. This suggests that the simulations 
contain the necessary physical ingredients to explain the 
observations, but the simulations do not by themselves 
reveal how these ingredients fit together to produce the 
observed result. Moreover, like the observed empirical 
rules, the simulations are based on a very limited range 
of galaxy properties, and in the absence of a model 
we can use to understand the origin of the simulation 
results, it is unclear how to extrapolate. Extending the 
simulations to cover the full range of galaxy parameters 
in which we are interested would be prohibitively 
expensive in terms of both computational and human 
time. 

Our goal in this paper is to remedy this lack of theo- 
retical understanding by providing a first-principles the- 
oretical calculation of the molecular con tent of a galactic 
disk i n terms of direct observables. In iKrumholz et al.l 
()2008l hereafter Paper I), we lay the groundwork for this 
treatment by solving the idealized problem of determin- 
ing the location of the atomic to molecular transition 
in a uniform spherical gas cloud bathed in a uniform, 
isotropic dissociating radiation field. In this paper we ap- 
ply our idealized model to atomic-molecular complexes in 
galaxies as a way of elucidating the underlying physical 
processes and parameters that determine the molecular 
content. We refer readers to Paper I for a full description 
of our solution to the idealized problem, but here we re- 
peat a central point: for a spherical cloud exposed to an 
isotropic dissociating radiation field, if we approximate 
the transition from atomic to molecular as occurring in 
an infinitely thin shell separating gas that is fully molec- 
ular from gas of negligible molecular content, the fraction 



of a cloud's radius at which this transition occurs is solely 
a function of two dimensionless numbers: 

_ fdiss^TdCEQ 

nTZ ^ ' 

TTi^nadR. (2) 

Here /diss ~ 0.1 is the fraction of absorptions of a Lyman- 
Werner (LW) band photons that produce H2 dissociation 
rather than simply excitation and radiative decay to a 
bound state, is the dust absorption (not extinction) 
cross-section per hydrogen nucleus in the LW band, Eq 
is the free-space number density of LW photons (i.e. far 
outside our cloud) , n is the number density of hydrogen 
nuclei in the atomic shielding layer, TZ is the H2 formation 
rate coefficient on dust grain surfaces, and R is the cloud 
radius. 

The quantity tr is simply a measure of the size of the 
cloud. It is the dust optical depth that a cloud would 
have if its density throughout were equal to its density 
in the atomic region. We may think of x as a dimension- 
less measure of the intensity of the dissociating radiation; 
formally it is equal to the ratio of the rate at which LW 
photons are absorbed by dust grains to the rate at which 
they are absorbed by hydrogen molecules in a parcel of 
predominantly atomic gas in dissociation equilibrium in 
free-space. This is a measure of the strength of the ra- 
diation field because in strong radiation fields the gas 
contains very few molecules, so most LW photons are ab- 
sorbed by dust and x is large. In weak radiation fields the 
gas contains more molecules, which due to their large res- 
onant cross-section dominate the absorption rate, mak- 
ing X small. 

Over the remainder of this paper we apply the model 
of Paper I to atomic-molecular complexes in galaxies. In 
§[2] we begin by considering giant clouds, which we may 
approximate as slabs, and in §[3] we extend our treatment 
to clouds of finite size. In § [4] we compare the model 
predictions of the previous two sections to observations 
of atomic and molecular gas. Finally in §[5] we summarize 
and discuss conclusions. 

2. THE ATOMIC ENVELOPES OF GIANT CLOUDS 

In this section we specialize to the case of giant clouds, 
which we define as those fo r which e'^^ ^ 1. For these 
clouds we show in iPaper II that the dust optical depth 
from the cloud surface to the atomic-molecuar transition, 
thi, is a function of x alone. We therefore begin our 
analysis with an estimate of x- 

2.1. The Normalized Radiation Field 

Of the quantities that enter into the normalized ra- 
diation field X, /diss is the most certain, because it is 
only a very weak function of the spectrum of the dis- 
sociating radiation. We therefore take it to have a 
constant value /diss ^ 0-1 independent of environmen t 
praine fc Bertoldilll996l : lBrowning erani2003l: Faper J ). 
Similarly, ad and TZ are functions of the properties of 
dust grains. These are both measures of the total sur- 
face area of dust grains mixed in the atomic shielding 
envelope around a molecular cloud; the former measures 
the area available for absorbing photons, while the lat- 
ter measures the area available for catalyzing H2 forma- 
tion. There will of course be an additional dependence 
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of these quantities on the optical and chemical proper- 
ties of grains, but these effects likely provide only a small 
fraction of the total variation. To first order, therefore, 
we expect the ratio a^/Ti- to vary little with galactic en- 
vironment, and we can simply adopt the value from the 
solar neighborhood. This is 
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where ad-21 — fd/lO^ cm^, TZ-ie.5 — TZ/IO 
cm^ s~^, and our best estimate s for the solar neighbor- 
hood give ua.-i\ ~ 7?-- 16.5 ~ 1 (jPraine fc Bertoldil[l996l : 
[VMfire ct al.' 2008'). 

Unfortunately, n and i?Q are considerably harder to 
determine, since we cannot easily make direct measure- 
ments of the atomic density around a molecular cloud or 
the dissociating radiation field to which it is subjected, 
particularly for clouds in extragalactic space. (It is pos- 
sible to determine these quantities for PDRs being pro- 
duc ed by individual sta r clusters - see lSmith et al.luOOOl 
and lHeiner et al.ir2008airbl - but these methods are gener- 
ally not able to determine mean radiation fields around 
giant clouds.) However, we can still gain considerable in- 
sight into the ratio / n that enters into x if we realize 
that n is not free to assume any value for a given Eq. 
The atomic gas in a galaxy generally comprises regions 
of both cold and warm gas (cold neutral medium and 
warm neutral medium, or CNM and WNM, respectively) 
in approximate pressure b alance (e.g. McKee_& Os trikcr 
119771 IWolfire et al.l[20M l. Molecular clouds form in re- 
gions where the gas is primarily cold. This is because 
the the effective opacity to LW photons provided by the 
small population of molecules found in a given element of 
predominantly atomic gas varies as n^, so the cold phase, 
due to its higher density, is far more effective at shielding 
from LW photons than the warm phase. Thus the n we 
are concerned with is not the mean density of a galaxy's 
atomic ISM, it is the density in the cold phase only. In 
the presumably dense gas in the vicinity of a molecular 
cloud, most of the mass is likely in the cold rather than 
the warm phase in any event. 

We can estimate the CNM density by using the con- 
dition of pressure balance between the cold and warm 
phases. 1 Wolfire et a l. (2003) show that, for a given am- 
bient FUV ra diation intensity Go (given in units of the 
iHabingl Il968l field, corresponding to a number density 
Eq 4.4 X 10~^ LW photons cm^^), ionization rate 
from EUV radiation and x-rays Ct, abundance of dust 
and polycyclic aromatic hydrocarbons Zd, and gas phase 
metal abundance Zg, the minimum number density nmin 
at which CNM can exist in pressure balance with WNM 
is well-approximated by 



l + 3.1(G[,zyCt'r^65 



(4) 



where the primes denote quantities nor malized to their 
values in the solar neig hborhood. .Wolfire et all (|2003l ) 
obtain this expression by constructing a temperature- 
density relation, determined by balancing the rate of 
grain photoelectric heating against cooling by the fine 
structure lines of Cii and Oi. Once they have con- 
structed the T ~ n curve, they identify the temperature 
at which the pressure is minimized. This is the warmest 



temperature at which the CNM can be in pressure bal- 
ance, and thus the corresponding density is the lowest 
possible CNM density. The primary uncertainty in this 
expression arises from the abundance, size distribution, 
and reaction properties of polycyclic aromatic hydrocar- 
bons (PAHs), but changes in PAH properties generally 
ch ange nmin only sA. t he factor of ^ 2 level (cf. Figure 8 
of lWolfire et"al]l2003l) . 

In a galaxy where young stars provide the dominant 
source of radiation and the IMF is constant, the FUV 
heating rate and the EUV/x-ray ionization rate are likely 
to be proportional to the star formation rate, and there- 
fore to each other. We therefore assume that Cj = Gq. 
Furthermore, if the physics of dust formation does not 
vary strongly from galaxy to galaxy, then the dust and 
gas phase metal abundances are likely proportional to 
the total metallicity Z, so we adopt Z'^ — Z'g — Z' . 
With these approximations the minimum CNM density 
becomes 
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cm 



(5) 



We caution at this point that both the assumptions that 
= Gq and Z'^ — Z'g — Z' are unlikely to hold in el- 
liptical galaxies, where young stars are not the dominant 
sources of EUV or x-ray radiation, and where the amount 
of dust per unit metallicity is known to be different than 
in spirals. Thus, equation is unlikely to hold in ellip- 
ticals. 

The CNM can exist in pressure balance at densities 
higher than nmin, so we take the typical CNM density to 
be 



'^CNM = (^CNM'^min- 



(6) 



We adopt (/)cnm ~ 3 as our fiducial value, which 
gives a CNM density of ncNM — 22 cm^'^ and (using 
IWolfire et al.l 's T — n relation, equation [TBI of this paper) 
a temperature of Tcnm = 105 K, consistent with obser- 
vations of typical CNM properties in the solar neighbor- 
hood. (Near a CMC, we expect Gg ~ 10 rather than 
Gq = 1 , due to the proximity of sites of star formation - 
Wolfire, HoUenbach, & McKee 2008, in preparation - but 
this does not affect our results, since we only care about 
the ratio G'q/ti.) In practice (/)cnm cannot be much larger 
than this, because pressure balance between the CNM 
and WNM is possible only over a limited range of CNM 
densities. If the CNM densities exceeds nmin by more 
than a factor of ~ 10 the CNM and WNM again cannot 
be in pressure balance because it is impossible for the 
warm phase to have a high enough pressure. 

Using equation Q for u^/TZ and equation ([6]) for n, 
and noting that the LW photon number densit y in the 
solar neighbor hood is roughly 7.5 x 10"'' cm"^ (jDraind 
[19781 IPaper If ), we find a total estimate for the dimen- 
sionless radiation field strength 



X = 2.3 



gd,-2i \ 1 + 3.1Z^° 

7^-16.5/ 0CNM 



(7) 



Note that all explicit dependence on the dust proper- 
ties, the radiation field, and the atomic gas density have 
cancelled out of this expression. 

Dependence on the dust properties has dropped out for 
the simple reason explained above: Ud and TZ are both 
measures of the dust surface area, so their ratio is nearly 
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constant. We can understand the somewhat more subtle 
reason that dependence on radiation field and the atomic 
gas density cancel by examining the physics behind ex- 
pression (|4]). As noted already, the minimum possible 
density in the cold atomic phase of the ISM corresponds 
to the density and temperature at which the pressure 
reaches a local minimum. Because the dependence of 
the cooling rate on gas temperature is determined al- 
most entirely by the quantum mechanical constants and 
element abundances that determine the shapes of the Cii 
and Oi cooling curves, and the photoelectric heating rate 
is essentially independent of temperature, the tempera- 
ture at which this minimum pressure occurs is nearly 
fixed at « 240 K, and does not depend on the back- 
gr ound rad iation field (c.f. equation 34 and Appendix C 
of IWolfir~e t al. 2003). Thus the density minimum will 
simply be the density at which the temperature reaches 
« 240 K. Since the heating rate varies as uEq and the 
cooling rate as n^, it immediately follows that the density 
at which a fixed temperature is reached varies as n cx -Eq ■ 
This explains why u/Eq is nearly constant in the CNM. 
There is only a weak dependence on metallicity, which 
arises because the heating rate depends on the charge 
state of polycyclic aromatic hydrocarbons (PAHs), and 
this in turn depends weakly on metallicity. 

Before moving on, we should note that we have ne- 
glected the role of internal radiation in determining 
where a cloud changes from atomic to molecular. This 
is justified because most stars that contribute significant 
amounts of dissociating radiation are born in molecular 
clouds, but they do not stay internal to those clouds for 
very long. Most dissociating photons come from massive 
stars born in clusters that burrow their way out of their 
parent molecular clouds via their winds and Hii regions 
in only a few Myr. Thus, most of the dissociating radia- 
tion to which a molecular cloud is subjected is delivered 
externally rather than internally, even if it comes from 
stars born in that cloud. 

This does, however, raise another cautionary point. 
We have computed the CNM density based on an im- 
plicit assumption of pressure balance, and we must con- 
sider under what circumstances pressure balance might 
not hold. One situation in which gas might not reach 
pressure balance is if it is subjected to hydrodynamic 
perturbations such as supernova shocks that create rapid 
and substantial changes in pressure, pushing gas into the 
unstable regime of density and pressure. Such gas is 
subject to an instability in which pockets of stable CNM 
condense within it, leaving behind a lower density ambi- 
ent medium that expands to become stable WNM (e.g. 
I Audit fc Hennebelldl2005l ). Thus the typical gas density 
will be significantly different than the value ticnm that 
we have estimated only if the time between successive 
shocks that drive gas into instability is small compared 
to the time required for this instability to operate. This 
is of o rder the cooling time scale, which IWolfire et al.l 
(|2003f l estimate to be 

1.2 / ^ X -0.8 



^cool ~ ^ 
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nT 



Myr, (8) 



104 Ky VSOOOKcm"^ 
where n and T are the gas number density and temper- 
ature. For typical CNM conditions near the Solar circle 
this is ~ 40 kyr, while for typical WNM conditions it is 
slightly under 10 Myr. Higher values of Gq, as are ex- 



pected near GMCs, reduce these to a few kyr and a few 
Myr, respectively. Thus perturbations that produce ve- 
locities 5 km s~^, i.e. not fast enough to induce shocks 
in the WNM, are essentially ineffective at driving the gas 
out of equilibrium, since any CNM gas they disturb will 
re-equilibrate very quickly. Stronger perturbations such 
as supernova blast waves that drive WNM unstable can 
keep substantial amounts of gas out of pressure balance 
only if they recur on time scales of a few Myr or less. 
For both supernova blast waves and the shocks induced 
by the ubiquitous turbulence in the atomic ISM the re- 
currence time is ~ 10 Myr, lar ger (although not h ugely 
so) than the equilibration time (jWolfire et al. I [20031 ) . We 
can therefore conclude that the typical atomic envelope 
around a molecular cloud is likely to be close to pressure 
balance between CNM and WNM. Some envelopes will 
have been subjected to a strong shock in the last few 
Myr, and these may have CNM densities substantially 
different than rtcNM, but they will be in the minority. 

Alternately, one could consider a galaxy in which the 
average galactic environment is so extreme that no two- 
phase equilibrium is possible, i.e. the pressure is so low 
that only WNM is stable, o r the pre ssure is so high that 
only CNM is stable. Wol fire et all and Schayc (200|) 
both find that a pressure so low that only WNM exists is 
consistent with vertical hydrostatic balance only in the 
very diffuse outer p arts of galactic disks. For the Milky 
Wav. IWolfire et all find that a CNM phase can exist ev- 
erywhere the mean ISM density is ^ 0.2 cm~'^, which is 
everywhere in the Milk y Way inside ~ 15 kpc in Galacto- 
centric radius: ISchavd estimates that hydrostatic balance 
requires the existence of a CNM phase any time the lo- 
cal gas surface density exceeds 3 — 10 Mq pc~^, which 
is true over a similar region. Thus, we can conclude that 
a pure WNM is unlikely to exist anywhere except in the 
far outer regions of galactic disks; in such regions, equa- 
tion ([7]) underestimates the dimensionless radiation field 
strength, and thus we will overestimate the molecular 
fraction. The converse possibility is a galaxy in which 
the pressure is so high that no WNM is present, only 
CNM; however, IWolfire et cdl find that that WNM can 
be present any time the mean ISM density ^ 70 cm~^, 
unless the CNM is entirely confined by a surrounding hot 
ionized medium. The case where no WNM exists because 
n 70 cm~'^ corresponds to an entire galaxy whose mean 
density matches that of a typical molecular cloud in the 
Milky Way, and such conditions are only found in star- 
burst systems where the molecular fraction is essentially 
unity. In such cases our models will provide only an up- 
per limit to the amount of Hi present, but even then 
our model may apply if the atomic shielding layer oc- 
curs far from the galactic midplane where the density is 
lower. The latter possibility, of an ISM consisting solely 
of cold atomic gas and hot ionized gas, appears not be 
be realized in nature. Thus, in summary, we expect our 
two-phase model to apply everywhere in galaxies except 
in their far outer parts, where the density is so low that 
no cold phase exists, and at the midplanes of starburst 
galaxies, where the density is so high that no warm phase 
can exist. 



2.2. The Shielding Column 
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The normalized radiation field x is the primary fac- 
tor controlling the size of the atomic gas column that 
is needed to shield a molecular cloud against dissocia- 
tion, and thus in determining fraction of the gas in a 
galaxy is molecular. This quantity in turn determines 
what fraction of a galaxy's ISM is molecular, and there- 
fore available for star formation, and what fraction is 
atomic, and thus without star formation. The invari- 
ance of x across galactic environments has important 
consequences for atomic-molecular complexes. First, x 
measures the relative importance of dust shielding and 
H2 self-shielding, so our results that x ~ 1 across galax- 
ies implies that dust- and self-shielding contribute nearly 
equally in essentially all galactic environments. We do 
not expect to find clouds where either dust- or molecular- 
shielding completely dominate except near strong sources 
of dissociating radiation where the atomic ISM is not in 
pressure balance or in regions of such low or high density 
that the atomic ISM does not have two phases. 

Second, the dust optical depth through the Hi shielding 
layer, thi, is primarily a function of x; the dependence on 
tr arises from geometric effects, and is greater than order 
unity only if the molecular region inside a cloud is a small 
fraction of its size. (Strictly speaking this is the optical 
depth of the CNM, not of all the atomic gas. As noted 
above most of the atomic gas around a molecular cloud is 
probably CNM, but there could be significant amounts of 
WNM along the line of sight that is not associated with 
that cloud, and which does not contribute to shielding 
it.) This implies that the dust optical depth through the 
atomic envelopes of molecular clouds should be roughly 
constant across galactic environments, at least as long as 
this optical depth is not close to that of the entire atomic- 
molecular complex. The only significant variation will be 
a weak increase in thi with metallicity. Correspondingly, 
the total Hi gas column will decrease with metallicity 
to a power less than unity, since thi increases slightly 
with metallicity, but the column of Hi required to achieve 
this dust optical depth decreases with metallicity. To 
illustrate for these effects we solve for thi as a fimction 
of Z in the limit tr — > 00 using the formalism of Paper I, 
and plot the result in Figure [TJ 

As the plot shows, four our fiducial model ^cnm — 3, 
we predict that the Hi layer around a molecular cloud in 
the Milky Way, Z' = 1, should have an absorption op- 
tical depth of Thi = 0.40 to LW photons, corresponding 
to an Hi column A^hi = 4.0 x lO^*' cm^^ (mass column 
density E = 4.5 Mq pc^^), assuming a dust absorption 
cross section per H nucleus of ~ 10~^^ cm~^. It is 
important to note that all of these values represent the 
absorption column on one side of a giant cloud. A 21-cm 
observation would detect the shielding column on both 
sides for a cloud exposed to the ISRF on both sides, so 
the detected Hi column would be double the values given 
in Figure [H As a shall see in §[3l the column is somewhat 
larger for a cloud of finite size. 

It is important to note that the shielding column we 
have calculated is somewhat different than the atomic- 
to-molecular tr ansition column den sity reported for the 
Milky Way by iSavaee et al.l (IToTTI logA^(g) = 20.7) 
and for the LMC and SMC bv iTumlinson et all ()2002l 
\ogN{H) > 21.3 and > 22, respectively). These values 
are the total column densities along pencil-beam lines of 
sight at which the fraction of the gas column in the form 
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Fig. 1. — The dust optical depth of the Hi shielding layer rjji 
{upper panel), and the corresponding Hi column A''hi {lower panel) 
for very large clouds, tjj^ 00. We show these results computed for 
"J^'CNM = 3 {thick lines) and for (f>CNM = 1 and 10 {upper and lower 
thin lines, respectively). The circles indicate our values for the 
Milky Way, Z' = 1, for our fiducial </)cnm = 3: thi = 0.40, A''hi = 
4.0 X 10^0. To compute Nm from Thi I w6 assume a dust absorption 
cross section per H nucleus in the LW band of crj = 10~^^ Z/Zq. 
We also show the visual extinction Ay corresponding to our thi 
and the mass column density S corresponding to A^hi- For the 
former we have assumed Av/rm = 0.48, following the models of 
IDraind (20033 b c) as explained in the text. To compute the latter 
we assume a mean particle mass per H nucleus of 2.34 X 10~^* g 
cm~^, corresponding to a standard cosmic mixture of H and He. 

of H2 reaches about 10% of the total. In cont rast, in 
the two-zone approximation we adopt in Paper I, we as- 
sume that the atomic-to-molecular transition is infinitely 
sharp, and under this approximation the shielding col- 
umn we report is the column at which the gas goes from 
fully atomic to fully molecular. Were the transition truly 
infinitely sharp as we have approximated it to be, the 
ratio of H2 to total column density would be zero at 
our computed shielding column. Comparing our theoret- 
ical shielding columns to the d etailed n umerical radiative 
transfer models we present in iPaper II shows that in re- 
ality, for conditions typical of the Milky Way, the ratio 
of H2 column to total column at our calculated transi- 
tion column A'hi is roughly 20%. Since this is a factor of 
2 larger than the 10% ratio used in the observationally- 
defined transition column, and in our simple model the 
H2 fraction increases linearly with total column density 
once we pass our predicted transition point, we expect 
our shielding column to be a facto r of ^ 2 larger than 
the values reported bv lSavage et al.l and ITumlinson et al.l 
We compare our model predictions to these data sets in 
more detail in § 14.21 

Using the extinction and absorption curves of iDrairi^ 
(|2003aibl H). the ratio of visua l extinc tion to 1000 A ab- 
sorption is Av/rm = 0.48 for iDraini 's Ry — 4.0 model, 
so the visual extinction corresponding to thi = 0.40 is 
Ay = 0.19. Adopting the Ry = 5.5 curve instead, 
appropriate for denser clouds, gives Ay = 0.28, while 
Ry — 3.1, for diffuse regions, gives Ay — 0.13. Our 
estimates for the LW dust optical depth and visual ex- 
tinction vary little with metallicity, changing by only a 
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factor of 2.7 for a metallicity ranging from IQ-'^Zq to 

The variation between the curves with 0cnm — 1, 3, 10 
show the full plausible range of variation in shielding col- 
umn arising from our uncertainty about the true den- 
sity in the atomic envelopes of molecular clouds. The 
0CNM = 1 and 10 curves are both within a factor of 2.6 
of the fiducial model, so this is an upper bound on our 
uncertainty. The actual error is likely to be smaller than 
this, since (/)cnm = 1 and 10 correspond to the extreme 
assumptions that the CNM assumes is minimum or max- 
imum possible equilibrium densities. 

We can obtain a quick approximation to the results 
shown in Figure [T] simply by noting that at solar metallic- 
ity our fiducial normalized radiation field is x = 3.1, and 
we show in Paper I that for a giant cloud with x < 4.1 
(corresponding to Z' < 2.5 for our fiducial parameters) 
the LW dust optical depth through the atomic shielding 
layer is 



THI = J, 



where 



= X 



2.5 



X 



2.5 + xe 



(9) 



(10) 



The dust-adjusted radiation field ip is a function only of 
metallicity; for our fiducial parameters 4>G'nm = 3 and 
o'd,-2i/'^-i6.5 = and Milky Way metalhcity, Z' = 
1, we obtain ip = 1.6. Moreover, the dependence on 
metallicity is weak: at Z' — 1/10, ip — 1.0, while at Z' = 
1/100, V' = 0.77. Because ^ depends only on metalhcity, 
we can also express the characteristic Hi shielding column 
on one side of a giant cloud solely as a function of Z': 



= 4.5 Af© pc 



fO,-21-Z'' 



(11) 

(12) 



where /xh = 2.34 x 10 g is the mean mass per hydrogen 
nucleus and CT0.-21 is the dust absorption cross-section at 
Milky Way metallicity (log Z' ~ 0) in units of 10^^^ cm^. 
The function f{Z', cf)^) is given by 



/(Z',(/)x) = 0.54 (0.32 + Z 



,0.385) ^-1 



where 



1.05(/)x-H0.42-hZ'0-385 
O.390X + 0.42 + 

0'd,-21 



, (13) 



(14) 



/(1, 1) = 1, and for our fiducial parameters 0^ = 1. The 
numerical factors that appear in f{Z',(j)^) are derived 
simply by substituting equation ([7]) for x into equation 
(dni) and thence into equation (|lip . These equations, 
and therefore the numerical values in the function /, de- 
pend solely on microphysical constants that describe the 
properties of molecular hydrogen and the chemistry of its 
formation on grain surfaces (which set /diss and cTd/TZ) 
and the shapes of the Cii and Oi cooling curves (which 
set the ratio E^/n). We have therefore calculated the 
shielding column to good approximation solely in terms 
of microphysical constants. 



3. THE ATOMIC ENVELOPES OF FINITE CLOUDS 

3.1. Formulation of the Problem 

To account for the fact that clouds have finite sizes 
and column densities, and that these can be quite small 
in dwarf galaxies or other low-pressure environments, 
we must examine the second dimensionless number that 
characterizes H2 formation and shielding: tr = na^R. 
Consider a cloud of known, fixed column density Scomp- 
If atomic-molecular complexes were of uniform density 
then we could find tr simply by multiplying Scomp 
by the dust cross section per unit mass ad/fJ-n, where 
/iH ~ 2.34 X 10"^'' g is the mean mass per H nucleus. 
However, the atomic region is warmer and has a lower 
mean mass per particle than the molecular one, and thus 
has a correspondingly lower density. This reduces the 
dust optical depth through it. Since it is the density and 
dust optical depth through the atomic shielding layer 
that matters, we must estimate tr, using the value of 
n appropriate for the atomic gas rather than the mean 
density in the complex. In other words, the quantity we 
want is 

Tr = ncNMO'd-R, (15) 

and we define 

(/'mol = (16) 
ncNM 

as the ratio of densities. Here nmoi and ncNM are the 
number densities of hydrogen nuclei in the molecular and 
CNM phases of the ISM, respectively. 

In the Milky Way, typical molecular cloud densities 
are n^oi « 100 cm"^ (McKee & Ostrikcr 2001), while 
observations of the giant Hi clouds around these molec- 
ular regions find typical dens ities ncNM ~ 10 cm"'^ 
(|Elmegreen fc ElmegreenI [19871 ). suggesting cjj^oi « 10. 
We do not expect this ratio to vary strongly between 
galaxies, so we should generally find 0moi ~ 10. 

We can can make an independent argument for (pmoi ~ 
10 by considering thermal pressure balance across the 
atomic-molecular interface. This argument only applies 
to gas near an atomic-molecular transition surface, which 
may or may not include the bulk of the gas in a cloud, 
but it does provide an estimate for the density ratio near 
the interface. Pressure balance requires that 

0mol = 1-8^- , (17) 

mol 

where Tmoi and Tcnm are the temperatures in the molec- 
ular and cold neutral atomic media, respectively, and the 
factor of 1.8 accounts for the difference in mean number 
of particles per H nucleus in the two phases. Across a 
very wide range of galactic environments the tempera- 
ture in the molecular phase of the ISM is T^oi ~ 10 — 20 
K, as a result of the balance between grain photoelectric 
heating and CO cooling. Since are interested in gas at 
the edge of the molecular region, we adopt Tmoi = 20 K 
as typical. (Temperatures are somewhat higher in star- 
burst systems, but in these galaxies the molecular frac- 
tion is essentially un ity in any event.) Using the model of 
IWolfire et al.l (|2003l ) for the atomic medium, and making 
the approximations Z'^ = Z'^ = Z' and (,[ = Gq as in 
§ 12. 1[ the density-temperature relation in the atomic gas 
is 

20G^rcNM.2e^-^/^'^-"'- 

J^CNM ~ 1 /, , (tSJ 



1 -I- 2.6(T( 



1/2 ;7/^0.365 
CNM, 2^ I 
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where Tcnm,2 = Tcnm/IIOO K). Combining this with 
equations (O and ([6|) enables us to write an imphcit 
equation for the CNM temperature in terms of 0cnm 
and Z': 



1.5/TcNM,2 



1 + 2.6(r^/fM.2^') 



^0.365 



0c 



31 



NM 



l + 3.1Z'0-365- 



(19) 



Substituting the solution to this equation in equation 
pT|) immediately gives us (/>moh the ratio of the number 
densities of H nuclei in the molecular and CNM gas. For 
our fiducial (/)cnm = 3, we find 0rnoi — 9.6 at Z' = 1, 
varying by only a few percent for metallicities in the 
range Z' = 10~^ — 10^. Given the encouraging agreement 
between this and the value (/>moi ~ 10 we find observa- 
tionally, we adopt the value of ^moi given by equations 
(fT7|) and (HH]) as our standard value for the remainder of 
this work. 

A spherical cloud that consists of a molecular core of 
number density rimoi and an outer atomic envelope of 
number density ncNM has a mean column density 



^comp — gMH^CNM^ [1 



fmol 



1)4. 



(20) 



where Xfja is the fraction of the cloud radius at which it 
transitions from molecular to atomic, i.e. = 1 cor- 
responds to a cloud that is molecular throughout and 
a^H2 = to one that is atomic throughout (see Paper 
I). It is convenient to rewrite this in terms of an optical 
depth 



3 

^ 4 



where E 



comp,0 



(21) 
(22) 

p/(1Mq pc~^) and the arrow in 



0.067 Z'Ecomp.o, 



the second step indicates that we have used our fiducial 

ad = 10-~^ ^ ^ ^ 



"^^Z' cm^. Equations ([20]) and imply 



= TR [1 



fmol 



-1)4.] 



(23) 



Note that neither tr nor Tc is the true center-to-edge dust 
optical depth of the complex; tr is the optical depth the 
complex would have if its density were ^cnm throughout, 
and Tc is the optical depth it would have if its atomic and 
molecular gas were mixed uniformly rather than spatially 
segregated. 

We are now in a position to compute the shielding 
column and the atomic and molecular fractions in fi- 
nite clouds. If we consider a complex of a given column 
density Scomp and metallicity Z', and we take the dust 
opacity to be given by its Milky Way value adjusted for 
metallicity, = CTd,MW^' with (Td,MW — 10~^^ cm~^, 
then equations ([21]) and ([23]) give one constraint on the 
unknowns tr and xr. from pressure balance between 
the atomic and molecular ph ases. Dis sociation- formation 
equilibrium, as comp uted inlPaper \ gives a second con- 



straint. We show in iPaper II how to compute the value 
of XRa for a given tr and x- its value is given implicitly 
by the solution to equations (33) and (37), or (43) and 
(44), of that paper. Since we have already computed x in 
terms of the metallicity (equation [7]), a choice of Scomp 
and Z' fully determine the two unknowns tr and xr. ■ 
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Fig. 2. — Visual extinction Ay (top panel), Hi column den- 
sity Shi {middle panel), and H2 fraction Mn^/M (bottom panel) 
in finite cloud complexes, as a function of complex mass column 



density Ec 



(or number column density A'^). In each panel the 



curves shown are for metallicities Z' running from 10-2 - 10\ as 
indicated. The hatched regions centered on the Z' = curve indi- 
cate the range of models with ipCNM = 1 ~ lOi while all other curves 
are for our fiducial value </>cnm = 3. In the top panel, the asterisks 
on the curves indicate the value of Scomp for which the molecular 
fraction falls to 1/2. The curves to the left of this point should be 
regarded as lower limits on the visual extinction to the molecular 
region. In the middle panel the two parallel dashed lines indicate 
ratios of Mjji/M = 1/2 and 1. Since values of Afjji below 1/2 are 
lower limits, once our models curves cross the lower of these lines, 
they true solution could be anywhere between them. Similarly, the 
dashed line in the bottom panel corresponds to Mn^/M = 1/2, 
and when the model curves fall below this they should be treated 
as upper limits. 

3.2. Numerical Solution 

We can either solve this system of nonlinear algebraic 
equations numerically, or approximate the solution ana- 
lytically. We first show the results of a numerical calcu- 
lation for a variety of values of Scomp and Z' in Figure 
[3 Rather than giving tr and xr. directly, which are 
not particularly interesting because we cannot measure 
them directly, we plot three derived quantities of interest 
which in at least some circumstances we can observe: the 
LW optical depth and visual extinction from the cloud 
surface to the atomic-molecular transition surface along 
a radial path, the Hi column averaged over the entire 
cloud, and the total II2 mass fraction over the entire 
cloud. 

It is important to point out that the value of Ay we 
report is measured differently than the Hi column den- 
sity Sri, or than the column densities we will use in 
§ [4| The column density is measured by averaging the 
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mass per unit area over the entire complex, while Ay 
is measured along a single pencil beam from the surface 
of the cloud to the atomic-molecular transition along a 
radial trajectory. The former quantity is more analogous 
to what is measured in an observation using a telescope 
beam that only marginally resolves or does not resolve 
a complex, while the latter is more closely analogous to 
a measurement of the extinction of a background point 
source through a cloud. Wc also caution that, for reasons 
we discuss in § 13.31 our predictions are only accurate for 
molecular mass fractions > 1/2. (This is in the worst 
case of very low metallicity and intermediate <?!>cnm ; our 
accuracy range expands as metallicity increases toward 
solar and as 4>cnm gets smaller or larger than 3.) Below 
this limit our calculations yield only upper limits on the 
molecular fraction, not firm predictions. This confidence 
limit is shown in the Figure [H 

The plots immediately yield a number of interesting 
results. First, our prediction of nearly constant Ay 
through the atomic shielding layers around molecular 
clouds continues to hold whenever there is a significant 
molecular fraction, even for finite clouds. Our prediction 
of a characteristic Ay ~ 0.2 through atomic shielding en- 
velopes of molecular clouds therefore continues to apply. 

We also find that there is a saturation in the Hi column 
density at roughly 6 Mq pc~^ for solar metallicity, which 
rises by a factor of a somewhat less then ten for every 
decade by which the metallicity declines. The Hi column 
saturates simply because once Ecomp is large enough, the 
cloud column densities become so large that they are ef- 
fectively in the infinite cloud limit. At this point the 
shielding column is geometry- independent, and is deter- 
mined solely by the normalized strength of the radiation 
field, Xj a- value that does not vary much from galaxy to 
galaxy. Once Scomp is sufficiently large to put a complex 
in the large cloud limit, adding additional mass simply 
increases the size of the shielded molecular layer, so the 
H2 fraction just rises smoothly. The saturation value of 
6 Mq pc~^ at solar metallicity is set by a combination of 
the fundamental constants describing H2 formation and 
dissociation, the shape of the Cii and Oi cooling curves 
(which determine the CNM density and temperature), 
and the properties of interstellar dust grains, which set 
the ratio a^/TZ. 

3.3. Geometric Uncertainties for Finite Clouds 

In § 4.7 of lPaper^ we show that our method for deter- 
mining the molecular abundance in finite clouds suffers 
from a systematic uncertainty arising from our imperfect 
knowledge of the opacity along rays that pass through the 
atomic envelope of a cloud. For our fiducial model we 
take the opacity along these rays due to molecules mixed 
into the atomic gas to be set by the value of the disso- 
ciation radiation field at the surface of the zone where 
molecules dominate the opacity. Our results depend on 
this approximation very little except at low molecular 



volume fraction. 



< 



0.2; in that case the uncertainty 



about this approximation means that our model enables 
us to predict only an upper limit on the molecular frac- 
tion, not an exact value. 

In this paper we are concerned with molecular mass 
rather than volume fractions, so we must quanti fy that 
uncertainty. To do so, we proceed as in § 4.7 of .Paper J : 
we adopt the opposite assumption, that opacity along 
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Fig. 3. — Ratio of minimum to maximum predicted H2 mass 
fraction versus maximum predicted H2 mass fraction, for a variety 
of values of Z' and (/"CNMi ^ indicated. Curves for logZ' < — 1 
are not shown because they are indistinguishable from those for 
logZ' = -1. 

rays passing through the region of the cloud where 
molecules dominate the opacity but still constitute a 
small fraction of all H nuclei is infinite. We then re- 
peat the calculations of § [3] following this assumption: 
for a given Scomp and Z', rather than solve the system 
of equations formed by equation P5)) of thi s paper and 
equations (33) and (37) or (43) and (44) of iPaperl . we 
instead solve equation (l23l) together with equations (69) 
and (70) or (71) of Paper I. Doing so gives a lower bound 
on the molecular content for a given Ecomp and Z' . By 
comparing the results in this case to our fiducial calcu- 
lations as presented in § [3l we obtain an estimate of the 
uncertainty of our results. 

In Figure [3] we show the results of this exercise. On 
the y-axis we show the ratio of the H2 mass fraction pre- 
dicted using our maximum opacity assumption, which 
represents the minimum possible molecular content, di- 
vided by the value produced by our fiducial model, which 
represents the maximum. This gives an estimate of our 
uncertainty. The x-axis indicates the H2 mass fraction 
predicted using the fiducial assumption we make else- 
where in the paper. As the plot shows, the two calcula- 
tions differ most at low Z' and intermediate (/icnm- In 
this case the calculations differ by a factor of a few for 
H2 mass fractions below ~ 0.5. If we adopt a factor of 
3 as an accuracy goal, this means that for cases where 
we predict an H2 mass fraction below 0.5, and at low 
metallicity, our predictions should be taken only as up- 
per limits. For solar metallicity or higher our confidence 
range extends down to molecular mass fractions around 
0.4, and we attain upper limits below this. 

3.4. Analytic Approximation 

We can gain additional insight into the behavior of the 
solution by constructing an analytic approximation. The 
ratio of the molecular mass AfH2 to the total complex 
mass M is 



/h2 



^inol2;H2 



M 



fmol 



(24) 



We wish to obtain an approximation for this in terms of 
the known quantities ip (given in terms of metallicity by 
equations [7] and llOp and Tc (given in term s of com plex 
column density by equation [2T|) . We show in lPaper^ that 
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for 7/1^3 and molecular volume fractfons x"^^ -^0.15, a 
range in parameter space that includes most of our mod- 
els for realistic cloud parameters, the molecular volume 
is well-approximated by 



1 



(25) 



4tra 

where for convenience we have defined 

tra -Tu + aTp, (26) 

and a = 0.2 is a numerical parameter that is optimized 
for agreement between the approximate and numerical 
solutions. Substituting this approximation into the con- 
dition for pressure balance, equation (|23p . gives 



»mol 



4tra 



(1 - 0mol) 



(27) 



As in lPaper^ . our approach to obtaining an analytic so- 
lution is to perform a series expansion in a. We therefore 
define 



TcA = Tc 



1 + ^ 



(28) 



so that Tca/tc = tra/tr- Using this definition of Tca 
together with equation (|27|) for Tc implies that 



TcA = 0mol Tra 



(29) 



If we now use our approximation (j25p and rewrite the 
result in terms of TcA using equation (|29p . we obtain 



/h2 



1 



4rcA 



(30) 



We must now express t^a in terms of ■0, Tc, and a alone. 
Thus 



TcA = Tc 



■ a'0 



■ a-ip 



TR 
\TRA 



(31) 



(32) 



The second term on the RHS still involves the unknowns 
TcA and Tra , but because they are already multiplied by 
a we now need only determine them to zeroth order in 
a. Solviing equation (p9)l for tra gives 



Tra 

TcA 



1 



Omol 
1 



□1 



1 



1 - 



1 ^ 


3V; 


0moi ; 


4TcA 


1 ^ 




0moiy 


4Tc' 



(33) 



where in the second step we have dropped a term of order 
a to obtain an expression that is accurate to zeroth order 
in a. Substituting this into equation (j32p . and thence into 
equation pop , gives our final expression for the molecular 
mass fraction, accurate to first order in a: 



37/; 
4^ 



1 



4a4"pi 



mol 



4rc -t- 3(0n 



1)V>, 



(35) 



Comparison of this approximate expression with the nu- 
merical solution illustrated in Figure [2] shows that for our 
fiducial 0cnm = 3 and metallicities from Z' ~ 10^^ — 10, 



it is accurate to better than 30% whenever the approx- 
imation analytic solution gives > 0.25, but that it 
goes to zero too sharply at low molecular fraction. We 
can improve the approximation by forcing the H2 frac- 
tion to approach zero smoothly rather than sharply at 
low column density. Experimentation shows that the ex- 
pression 



47^ 

3^- 



1 



mol 



4Tc + 3(0,nol - 1)V' 



(36) 



matches the numerical result for /hi = Mhi/M for 
0CNM = 3 to better than 20% for all Z' < 10 regardless 
of the value of /hi. (However note that, as we show in 
§ 13.31 for /h2 ^1/2 our estimate of the molecular content 
is only an upper limit, and this is true of equation 1361 as 
well.) Using equation (|2T|) to replace 
substituting in our fiducial values 4>cnm 
and (Td = IQ^^^Z' cm^^, equation 



with Scomp, and 
= 5, a = 0.2, 
becomes 



rai 



^ , s \3 /125 + s 



-1/3 



where 



and Y.r 



^comp 



4, 



(37) 



(38) 



^comp.o = Scomp/(l Mq pc~ ). Notc that our result 
indicates that to good approximation the molecular con- 
tent of an atomic-molecular complex depends only on the 
combination of input parameters Z'I]comp,o/'0; the nu- 
merator Z'Ecomp.o is simply the dust column density of 
the complex up to a scaling factor, while the denomina- 
tor ijj is the dimensionless radiation field, which equations 
([7|) and (fTO| give solely as a function of metallicity. 

From (l37|) . it also immediately follows that the H2 to 
Hi ratio Rb_^ = /H2/./HI is 



/125 



V 96-l-s 



1/3 



1. 



(39) 



For i?H2 > I7 which is the regime for which our models 
apply with high confidence, we can use an even simpler 
expression 



i?H2 



0.08s = 0.08- 



-'comp 



(40) 



where is accurate to ~ 30%. 

Similarly, we show in lPape^ that the dust absorption 
optical depth through the atomic layer for a finite cloud 
is well-approximated by 



Tffl ~ -r 



1 



1 - (a'/4)(VVTR) 



(41) 



where a' = | — 4a = 0.7. If we treat a'/4 as a small 
parameter and perform a series expansion around it, then 
we need only approximate tr to zeroth order in a. We 
can do this simply by using equation (|34p with a = 0, 
which allows us to set tr = tra and Tc = Tca- Thus to 
zeroth order in a we have 



Tr — Tra 



1 

^mol 



1 - 



1 



»inol 



3lp 
4^ 



(42) 



10 



and to first order in a or a' we have 



1 



a'llxf) 



mol 



4re + 3(0„ 



(43) 



As with approximation ((35|) for the molecular mass frac- 
tion, this expression works well whenever the molecular 
fraction is not too low, and may be improved by forc- 
ing the optical depth to approach the total cloud optical 
depth smoothly when the column density becomes low. 
The expression 



'HI 



16 



1 



4t, 



(44) 



3((/),„oi - 1)V'. 

is accurate to better than 35% for all Z' < 10, and to 
better than 25% for Z' < I. 

It is also convenient to invert our analytic expressions 
to determine column density as a function of molecular 
content and metallicity. The term (125 -I- s)/(96 + s) in 
equation (|37|) is generally close to unity except at ex- 
tremely high column densities, so if we neglect second- 
order corrections to the difference between this term and 
unity, we can solve equation (j37|) for /hi to obtain 



■^comp 



HI^ 



1 



1/3 



8.7+(/Hf-l 
11+(/hi'-1) 



,1/3 



.1/3 



(45) 



This expression matches the numerical solution at the 
20% level for /hi < 0.75. Note that this result im- 
plies that the column density Scomp at which a given Hi 
fraction is reached depends on metallicity both explic- 
itly through the Z' term in the numerator, representing 
the effect of metallicity on dust content, and implicitly 
through ip feguations [71 and [TU)). representing the effect 
of metallicity on the ratio of radiation intensity to CNM 
density. Since "0 is an increasing function of metallicity, 
for a given molecular fraction Scomp has a weaker than 
linear dependence on metallicity. For example, evalu- 
ating with /hi = 0.5 at solar metallicity Z' = 1 
indicates that we expect the gas to be half molecular 
for complexes with Scomp — 27 M© pc~^. (The exact 
numerical solution is Scomp — 25.5 AIq pc~^.) At one- 
third solar metallicity, Z' = 1/3, half molecular content 
is reached at Scomp — 67 Mq pc^^ (using equation [45{ 
numerically Scomp = 55.3 AIq pc~^), somewhat less than 
a factor of 3 higher. 



4. COMPARISON TO OBSERVATIONS 

Our model makes strong predictions for the relative 
fractions of Hi and H2 as a function of total surface den- 
sity and metallicity, and in this section we compare these 
results to a variety of galactic and extragalactic observa- 
tions. 

4.1. Extragalactic Observations 
4.1.1. Data Sets 

We use three extragalactic data sets for comparison to 
our models. Two are recent surveys that have mapped 
nearby galaxies in 21 cm Hi and 2.6 mm C0(1 0) 
emission at overlapping positions, and therefore provide 
an ideal laboratory in which t o test o ur model. The first 
of these is the work of WBOtlBROl . and ,BR06„ who re- 
port Hi and H2 surface densities on a pixel-by-pixel basis 
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dicatcs that there is no gas-phase metallicity is reported in the 
literature. 

We take the metallicity relative to solar to be proportional 



to the O/H ratio, i.e. logZ' 
galaxy is in ,BROd sample, L : 



: [log{0/H) + 12] - 8.76. B 
: galaxy is in IL08I sample 



in 14 nearby galaxies (including the Milky Way). The 
H2 surface densities are inferred from CO observations 
taken as part of the BIM A SONG survey (| Regan et al.l 
I2OOII: iHelfer et al.]|2003l) . while t he Hi observations are 
from the VLA. lWB02l "and IBR06I supplement these data 
with stellar surface density measurements from 2MASS 
(Jarrett ct al. 2003), which together with the equations 
given in 1BRO6 can be used to derive a pressure in each 
pixel if one assumes that the gas in in hydrostatic bal- 
ance, that the stellar scale height greatly exceeds the gas 
scale height, and that the gas velocity dispersion has a 
known value. The galaxies in the sample are all molecule- 
rich spirals with metallicities within 0.5 dex of solar. 

Th e second extragalactic data set we use is compiled by 
iLOSl who a combine Hi measurements from the THINGS 
survey with CO data partly taken from BIMA SONG and 
partly from the ongoing HERACLES survey. The au- 
thors also include 2MASS stellar surface densities in their 
compilation, and give mean pressures. Unlike the .BR06I 
sample the data reported are averages over galactocen- 
tric rings rather than individual pixels, although point- 
by-point maps at sub-kpc resolution are in preparation 
(F. Walter, 2008, private communication). The sample 
includes 23 galaxies, of which roughly half are large spi- 
rals and roughly half are low-mass, Hi-dominated dwarfs. 
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The galaxies i n the d ata set partly overlap with those in 
the sample of IBROGI but extend over a wider range of 
metallicities and molecule fractions. 

We summarize the galaxies in the samples in Table 
[H We also report metallicities for each galaxy where 
these are available in the literature. We do not include 
NGC4414 in the analysis, because no gas phase metal- 
licity is available for it. In the comparison that follows, 
we neglect the presence of metallicity gradients within 
these galaxies, because gradients are only available for 
some of them. On top of this, we note that the metal- 
licities themselves are probably uncertain at levels from 
hundreths to tenths of a dex, depending on the galaxy 
and the analysis technique used, which adds additional 
scatter on top of that already introduced by our neglect 
of metallicity gradients. 

In both data sets the uncertainty in the Hi column 
densities is generally ~ 10%. In the CO data the formal 
uncertainties are generally ~ 20 — 30%, but the dom- 
inant uncertainty is probably systematic: the X factor 
used to convert observed CO luminosities into II2 col- 
umn densities. This i s uncertain at the factor of 2 
level (iBlitz et al.ll2007l), and almost certainly varies with 
metaUicitv fe.g. iBolatto et al.ll2008f) . Thus, ahhough for 
clarity we will suppress error bars in the plots that fol- 
low, recall that the molecular data is uncertain at the 
factor of ~ 2 level. 

The third data set we use, in § 14.1.41 is the 
Syitzer Survey of the Small Magellan ic Cloud (S^MC) 
(|Bolatto et al.ll2007HLe"rov et al.ll2007D . This survey dif- 
fers from the SONG and THINGS data sets in that those 
surveys infer the presence of H2 via CO emission, whereas 
S^MC measures molecular hydrogen using measurements 
of dust from Spitzer comb ine d with Hi measurements b y 
IStanimirovic et all ()1999[ ) andlStanimirov ic et aT1(|2004l ). 
The basic idea behind the technique is that one deter- 
mines the dust-to-gas ratio in a low-column density re- 
gion where molecules are thought to contribute negli- 
gibly to the total column. Then by comparing the Hi 
and dust column density maps, one can infer the pres- 
ence of H2 in pixels where the total dust column ex- 
ceeds what one would expect for the observed Hi col- 
umn and a fixed dust-to-gas ratio. The reason for using 
this technique is that, at the low me tallicity of the SMC 
(log[0/H] + 12 = 8.0, lDufouil[l98l i.e. 0.76 dex below 
Solar), CO may cease to be a reliable tracer of molecular 
gas. The SMC represents the lowest metallicity galaxy 
for whi ch we have H2 detections rather than upper lim- 
its; the IBROGI sample does not contain any galaxies with 
metallicities as low as the SMC, and no CO was detected 
in any of the lL08l galaxies with metallicities comparable 
to or lower than the SMC. Thus, the SMC represents a 
unique opportunity to test our models at very low metal- 
licity. 

Before comparing to these data sets, it is worth com- 
menting briefly on one additional extragalactic data set 
to which we will not compare our models: observations 
of H2 column densi ties along sightlines in the LMC and 
SMC using FUSE (|Tumlinson et al.l[200l) . We do not 
use this data set for comparison because it includes only 
sightlines with low column densities that are strongly 
dominated by atomic gas; the highest reported H2 frac- 
tion is below 10%. Sightlines with significantly higher 
molecular content than this absorb too much background 



starlight to allow FUSE to make a reliable measurement 
of the H2 column. The low column densities of the clouds 
that FUSE can observe in these galaxies place almost all 
of them into the regime where our theory yields only up- 
per limits. Those limits are generally consistent with the 
data, but the comparison is not particularly illuminating. 

4.1.2. Column Density and Metallicity Dependence 

In this section we compare our models predictions of 
the Hi and H2 mass fractions as a function of total 
column density and metallicity to our two extragalac- 
tic data sets. To perform this comparison, for conve- 
nience we first break the sample into four metallicity 
bins: -1.25 < logZ' < -0.5, -0.5 < \ogZ' < -0.25, 
-0.25 < \ogZ' < 0, and < \ogZ' < 0.5, where 
Z' = Z/Zq and we adopt log(0/H) + 1 2 = 8.76 as 
the v alue corresponding to solar metallicity (jCaffau et al.l 
[200I . These bins roughly evenly divide the data; we do 
not include the galaxies for which metallicities are not 
available in the literature. 

Next we must consider how finite spatial resolution will 
affect our comparison. Given the beam sizes in the ob- 
served data, which range from ~ 0.3 kpc for the nearest 
gala xies to a few kpc for t he most distant, each pixel 
(for Hoi) or ring (for 'L08') in the observed data set is 
likely to contain multiple atomic-molecular complexes. 
The atomic and molecular observations are convolved to 
the same resolution, so this does not bias the measure- 
ment of the atomic to molecular ratio, although it does 
mean that this ratio is measured over an averaging scale 
set by the beam size. Since individual atomic molecu- 
lar complexes presumably represent peaks of the galactic 
column density, however, the total (Hi -I- H2) observed 
gas surface density Eobs reported for each pixel or ring 
represents only a lower limit on Scomp, the total sur- 
face density of individual complexes. In our model the 
fraction of the cloud in molecular form is a strictly in- 
creasing function of Scomp (and of metallicity Z), while 
the atomic fraction is a strictly decreasing function. Be- 
cause we have only the observed column density Eobs 
available, not the column density of an individual com- 
plex Ecomp, we have no choice but to take Scomp = Eobs 
as our best estimate. The error in this approximation 
probably ranges from tens of percent for nearby galax- 
ies where the beam size is not much larger than the size 
of a complex up to an order of magnitude or more for 
ring-averages in distant galaxies. As a result, though, we 
expect to overestimate the atomic fraction and underes- 
timate the molecular fraction. Physically we may think 
of this as a clumping effect: the clumpier the gas is, the 
better able it is to shield itself against dissociating radi- 
ation. Since our observations smooth over scales larger 
than the characteristic gas clumping scale, we will miss 
this effect. 

There are also two effects which go in the other di- 
rection, however. First, as noted above, shielding comes 
primarily from cold atomic gas, not warm gas. Although 
most of the gas in the immediate vicinity of a single 
molecular region is probably cold, observations that av- 
erage over many molecular cloud complexes are likely to 
include a fair amount of WNM gas as well. Since we 
have only the total Hi column densities including both 
phases, we are prone to overestimate Ecomp and therefore 
the molecular fraction, because the WNM raises the Hi 
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column but does not provide much shielding. This effect 
is not important at moderate to high column densities, 
where the atomic gas does not dominate the total mass 
budget, but it could be significant at lower densities. 

Second, there may be gas along our line of sight 
through a galaxy that is not associated with an atomic- 
molecular complex along that line of sight. This effect 
will increase in severity as the galaxy comes close to edge- 
on, since this will increase the path length of our line of 
sight through the galaxy. As with WNM gas, this ex- 
tra material contributes to Sobs but not to the complex 
surface density Scomp, and it therefore leads us to over- 
estimate the molecular fraction. We can perform a very 
simple calculation to estimate the size of this effect. Con- 
sider a simple self-gravitating gas disk characterized by 
the standard vertical density profile n oc sech^(0.88z//i), 
where h is the half-height of the gas. At some point 
in this disk is an atomic-molecular complex, centered at 
the midplane. Since complexes are found preferentially 
at the midplane, unless the galaxy is very close to edge- 
on then we need not consider the possibility of our line of 
sight intersecting multiple independent complexes. Since 
the complex formed from a large-scale gravitational in- 
stability in the disk, its characteristic size is ~ h, and we 
therefore consider gas to be "associated" with the com- 
plex if it is within a distance h of it in the plane of the 
disk. Suppose this galaxy has an inclination i. Mak- 
ing the worst-case assumption, that there is no density 
enhancement due to the presence of the self-gravitating 
complex, we can then compute what the fraction of the 
gas we see along our line of sight is not associated with 
it. This is simply 



/^°^^^.sech^(0.88z//i) dz 
sech^{0.88z/h)dz 



= 1 - tanh(1.13coti). (46) 



This is less than 0.5 for all inclinations less than 64°. 
Only a handful of the galaxies in the SONG and THINGS 
surveys have inclinations larger than this, so even in the 
worst case scenario where complexes do not represent any 
enhancement of the gas density, we expect non-associated 
gas to produce an error smaller than a factor of ~ 2 
for the great majority of the galaxies to which we are 
comparing. 

Given the limitations imposed by finite resolution, we 
proceed as follows. Using the method described in § [3l 
we compute the molecular mass fraction 

/h2 (^comp; ^) ^ 



M 



(47) 



and 



as a function of complex column density T,^ 
metallicity Z. (For these and all subsequent predictions 
we use our fiducial value of 0cnm = 3.) We then generate 
predicted Hi and H2 column densities for each metallicity 
bin and each observed column density E via 



(48) 
(49) 



.predicted — (^obs, •^min)5jobs 
Shi, predicted = [1 ^ (^obs, ^inin)]Sobs 

where Zmin is the minimum metallicity for that bin. 
Since /h2 (^comp, ^) is a strictly increasing function of 
^comp and Z, and we know that Sobs < Scomp and 
Z^in < Z, we expect /h2 (^obs, ^min) to be a lower 
limit on the true molecular fraction. We therefore expect 
Sh2, predicted to bc a lowcr limit on the observed data, and 
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Fig. 4. — Hi column density Ehi versus total column density 
S for galaxies in metallicity bins —1.25 < logZ' < —0.5, —0.5 < 
logZ' < -0.25, -0.25 < logZ' < 0, and < logZ' < 0.5, as 
indicated. In each p anel we plot the v alues of Sjji and S from 
the samples of IBROa (squares) and ILOSl (diamonds) and show our 
model predictions of Shi as a function of S for logZ' = —1.25, 
—0.5, —0.25, and (lines, highest to lowest). The curve for the 
value of log Z' equal to the minimum log Z' for each bin, which 
should represent the upper envelope of the data, is shown as a 
solid line. The rest are shown as dotted lines. The four lines 
are the same in each panel. To maximize readability we omit the 
error bars on the data points. The parallel slanted dashed lines, 
Sfji = Soi3s/2 and Sfji = Sq^s, show the range of value of Shi 
for which our predictio ns s hould be treated as upper limits for the 
reasons discussed in § 13.31 Parts of our model curves above the 
Sfji = Sobs/2 line could be as high as the Sjji = S^j^g line. 

Shi, predicted to be an upper limit. The possible exception 
to this statement is at low column densities, where a sig- 
nificant fraction of the Hi column may be in the form of 
WNM. 

We plot the data against our theoretical prediction for 
the upper envelope of Shi in Figure 21 and we show 
the corresponding predicted lower envelopes for Sh, and 
i?H2 = Shj/Shi in Figures Eland El For the |BR06| data 
set, rather than plotting the tens of thousands of individ- 
ual pixels it contains, for each galaxy we show the data 
averaged over 20 logarithmically-spaced column density 
bins running from the minimum to the maximum value 
of Sobs reported for that galaxy. As the plots show, our 
model predictions for the upper envelope of the Hi sur- 
face density and the corresponding lower envelope of the 
H2 surface density as a function of total surface density 
and metallicity agree very well with the data. The data 
fill the space up to our predicted envelopes but for the 
most part do not cross them, even when the predicted Hi 

envelope becomes fiat for Sobs ^ 10 Mq pc^^. Moreover, 
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Fig. 5. — Same as Figure |4l except that we show H2 rather than 
Hi column densities. For clarity we do not show non-detections. 
Note that AfH2 is the column density of H nuclei in molecular form, 
which is twice the column density of H2 molecules. The slanted 
dashed line is Sjjj = Soi,s/2; as discussed in ^ 13.31 above this line 
our model curves may be taken as predictions, while below it they 
should be taken as upper limits. 



our model recovers not only the primary dependence of 
Shi on the total observed column density Sobs, but also 
the secondary dependence on metallicity. For example, 
rings in the lowest metallicity bin in the iLOSl data set 
reach total mean column densities of almost 20 Mq pc~^, 
but still show no detectable molecular component. On 
the other hand, in the highest metallicity bin the molec- 
ular fraction is close to 70% in rings with Sobs ~ 20 Mq 
pc~^. Our models reproduce this effect: at a metallicity 
of logZ' — —0.5, we predict that the gas will be 94% 
atomic even at a surface density of 20 Mq pc~^, whereas 
for log Z' = 0.5 we predict an atomic fraction of only 
30% at that column density, in agreement with the data. 

We caution that we can only predict upper limits on 
the molecular content in regions of parameter space when 
our predicted molecular fraction falls below ~ 1/2, for 
the reasons discussed in § 13.31 We have indicated the 
regions where our model predictions convert to upper 
limits in Figures 2] - [S] Alternately, one can express 
this uncertainty as giving a minimum column density at 
which we can predict a value rather than an upper limit 
for molecular content. For reference, at the metallicities 
of logZ' = -1.25, -0.5, -0.25, 0, and 0.5 which define 
the edges of our metallicity bins, the minimum column 
densities for which we can predict numerical values to 
better than factor of few confidence are 250, 58, 38, 25, 
and 11 Mq pc~^, respectively. 
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Fig. 6. — Same as Figures |4] and O except that we show i?H2 = 
Sjjj/Shi rather than column densities. For clarity we do not show 
non-detections. Th e hor izontal dashed line corresponds to J?H2 = 
1; as discussed in 13.31 above this line our model curves may be 
taken as predictions, while below it they should be taken as upper 
limits. 



Finally, we note that the future HERACLES / 
THINGS data set represents an opportunity to perform 
an even stronger test of our model. In Figure [3] a signif- 
icant fraction of the data points fall below our predicted 
upper limits, and correspondingly these points are above 
our lower limits in Figures [5] and [HI We hypothesize that 
these data points represent rings or pixels within which 
the gas is significantly clumped, so that the averaged 
column density Sobs seen in the observation is signifi- 
cantly lower than the column density at which most of 
the molecular gas in that beam or ring is found, and our 
calculation for a complex with Scomp — Sobs overesti- 
mates the Hi column. In reality these regions probably 
consist of patches of high column density where most of 
the molecules reside, embedded in a lower density ambi- 
ent medium that has a lower molecular fraction than we 
determine by averaging over large scales. If we could ob- 
serve these regions at higher resolution, in Figure [4] the 
high column, high molecular content points will lie to the 
right of and slightly above the low resolution points, since 
both the total and Hi column densities will be higher 
than for the lower resolution observation, but the in- 
crease in the total column will be larger than in the 
Hi column. Conversely, the low column, low molecule 
patches will lie to the left and slightly downward from 
the low resolution points, since both the total and Hi 
columns will decline, but the Hi by less, since the atomic 
fraction rises. These changes will bring the data points 
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Fi g. 7. — Same as Figure |4] except that we show only data 
from lBROel . and we have divided this data set into "near" galaxies, 
those for which the survey resolution is < 1 kpc, and "far" galaxies, 
those for which the resolution is > 1 kpc. 



closer to our model curves. Indeed, the existing data 
already hint that such an effect i s prese nt: the single 
beam-averaged observations of the |BI106| data set scat- 
ter away from our limit lines n otice ably less than the 
ring-averaged observations from iL08|. Sim ilarlv. we can 
divide the point-by-point data from IbROGI into a "near" 
sample, consisting of galaxies for which the resolution 
is smaller than 1 kpc, and a "far" sample, consisting of 
galaxies with larger resolutions. We plot a v ersion of 
Figure H] using only this divided data set from IBR06I in 
Figure [T] The comparison is quite noisy, but the data in 
that bin do at least seem consistent with the hypothesis 
that the near data fall closer to the model lines than the 
far data. 

The full HERACLES survey, currently underway, will 
report measurements of the molecular surface density for 
individual patches 0.5 kpc in size, generally smaller 
than the beam patches the BIMA / SONG survey. While 

this is still considerably larger than the ^0.1 kpc-size of 
a typical atomic-molecular complex in the Milky Way, at 
the higher spatial resolution of the full HERACLES data 
set beam-smearing effects should be reduced and the col- 
umn densities reported for each patch should be closer to 
the true column densities Scomp of the individual atomic- 
molecular complexes. We therefore predict that in the 
full HERACLES data set the Ehi — Sobs relation should 
be closer to our theoretical upper and lower limit curves 
than the lower resolution or azimuthally-averaged data 
shown in Figures |4] - [6l 



4.1.3. Pressure Dependence 

IBR04I and IBR06I find that the molecular fraction in a 
galaxy correlates with the mid-plane gas pressure as 

0.92±0.07 



i?H2 



Shi V3.5±0.6 X 104 K cm 
while Eol find 

Ry{. 



2.0 X 104 K cm"^ 



0.8 



(50) 
(51) 



for their sample. Although our model does not directly 
give a prediction for the dependence of the molecular 
fraction on pressure, and we argue that surface density 
and metallicity are the physical variables that directly 
control the molecular fraction, wc wish to check whether 
our model is consistent with the observed correlation. 

To avoid introducing any bias in performing this check, 
we must determine p ressure s fr om o bservable quantities 
in the same way that lBR06l andjL08ldo. following an ap- 
proximation introduced bv IBRO j .^ The BR04 approxi- 
mation treats the galaxy as an infinite thin disk of uni- 
form gas and stars in vertical hydrostatic balance. For 
such a disk the pressure is related to surface density by 



P = 2 



(52) 



where Eg is the total (Hi plus H2) gas surface density, 
is the stellar surface density, and Vg and are the verti- 
cal velocity dispersions of the gas and stellar components, 
respectively. The gas velocity dispersion is roughly con- 
stant, tig ~ 8 km s~^, for the galaxies in the IBR06I and 
iLOSl samples. The term E*/w* presents more difficulty, 
however. It varies by orders of magnitude from galaxy 
centers to edges, and between galaxies, so we cannot pick 
a single value for it. Since in our model the molecular 
fraction is a function only of metallicity and Eg, this 
means that we do not predict a single-valued relation- 
ship between Rn^ and P that we can directly compare 
to the empirical fits. 

However, we can still compare our model to the data in 
two ways. First, we can pick a range of values of E^/w^^ 
consistent with the range in the observed sample, and 
demonstrate that the resulting range of predictions for 
i?H2 covered by our model is consistent with t he observa- 
tional data. For this purpose we follow lBR06l in adopting 
a constant value fg = 8 km s~^ and using simplified ver- 
sion of equation ([5^ . which follows from assuming that 
both the surface density and scale height of the stars are 
much greater than those of the gas: 



P 

— = 272 cm" 
kn 



K 



Mq PC-^J \Mq pC-2 
0.5 



0.5 



km s ^ / V PC 



(53) 



^ It is important to note that the midplane pressure which the 
observers attempt to estimate is not the same as the pressure in 
the molecular gas or the CNM which we balanced in § 12.11 the 
total midplane pressure includes contributions from all phases of 
the ISM rather than just the CNM or molecular clouds, as well as 
contributions from magnetic fields, cosmic rays, and bulk motions. 
The pressure we use in §[2TT] includes thermal pressure in the CNM 
or molecular gas only. 



15 




Fig. 8. — Molecular to atomic ratio -Rh2 versus pressure P. We 
show the data sets of BR06 (squares) and XOSi (diamonds), and 
predictions for models with logp* = —2.5, —1.5, —0.5, 0.5, and 
1.5 in units of Mq pc~^ (highest to lowest lines, as indicated), 
and with metallicities logZ' = —0.5 and (dashed and solid lines, 
respectively). In the upper panel for clarity we show only a single 
galaxy from each data set (NGC4736 from ,BR06i and HOII from 
ILOSI) . while in the lower panel we show the data for all galaxies; 
the model curves in each panel are identical. For IBR06I the data 
we show are averaged ov er 0. 1 dex-wide bins in pressure, as in 
IBROQ' s Figure 3. For the ILOSI data, we only plot rings for which 
both molecular and atomic gas are detected. The dotted horizontal 
lines show -Rh2 = Ij below this line our predictions should be taken 
as upper limits only. 

where Eobs is the observed value of the total gas surface 
density Eg, h^, is the stellar scale height, and for con- 
venience we define — We consider values of 
logp* from —2.5 to 1.5 in units of Mf7^ pc~ ^, a range that 
covers almost all of the samples of IBR06I and|L08. For a 
given value of p*, equation ([55)1 gives a single- valued re- 
lationship between P and Eg, so for each P and a choice 
of metallicity we can use our model to make a predic- 
tion for Ehi and thus for Rn^ ■ As with our calculations 
in § 14.1.21 we expect to systematically underpredict R^^ 
because the observed column density Eobs is an under- 
estimate of the true column density of atomic-molecular 
complexes that are not resolved by the telescope beam. 

We show the results of this computation in Figure [H 
As expected, the model covers an area consistent with 
the observed data. First note th at the re is a sys tematic 
offset between the results from IBR06I and lL08l : this is 
likely because some of the structural parameters that are 
used in equation ([55]) are quite uncertain, particularly 
/i*, and the values adopted by the two surveys are not 
necessarily the same. We should keep this uncertainty in 
mind as we proceed, since it suggests an upper limit on 
the level of agreement we can expect. 

Nonetheless, the correlation between i?H2 and P shown 



-0.9 



IS some- 



in both observational surveys, i?H2 P 
what flatter than any of our model curves for a particular 
choice of p, , which all approach Rn^ oc P^ at large P. 
However this is to be expected: low pressures and gas col- 



umn densities are found preferentially in the outer parts 
of galaxies where p* is small, so at low P we expect to 
be closer to the higher model curves, which have low p,. 
Conversely, high pressures and values of Eg are found 
preferentially in regions with higher stellar densities p,, 
so high P values should be closer to the lower model 
curves, which have larger p*. This covariance between 
P and p* results in a flattening of the Rn^ oc P relation 
that we would predict if all galaxies had fixed p*. 

We can eliminate this covariance effect and produce a 
stronger test than that shown in Figure [8] by using our 
model to generate predictions for molecular content di- 
rectly, and plotting those against the inferred pres sure. 
To do this, for each galactocentric ring in the |L08| data 
set we take the observed total gas surface density Eobs 
and the metallicity given in Table [T] (which we treat as 
constant in each the gal axy) an d use our model to predict 
i?H2 ■ Similarly, for the IBR06I data set we take the ob- 
served gas surface density Eobs in each pixel and use our 
model to generate a prediction for R^^ . Since the pres- 
sure in each ring or pixel is known, by this procedure we 
generate a synthetic set of data points (P, Rn^ ) which we 
can compare to the observations, and from which we can 
generate a fit for i?H2 as a function of P. 

Before making the comparison, however, we must also 
model the effects of finite telescope sensitivity, which 
make it impossible to detect molecules below a certain 
minimum column density. Since, as Figure [5] shows, our 
model predicts that the slope of Rn^ versus P varies with 
P, this can affect the fit we generate from our synthetic 
data. We include this ef fect using a procedure nearly 
identical to that of lBR06l These authors first estimate 
the minimum value of i?H2 as a function of pressure that 
they could detect in each galaxy based on their telescope 
noise limits. For our synthetic data we adopt a minimum 
^H2.min = 0.05 for all ga laxies, slightly below the lowest 
average R^^ that lBR06l report for any of the galaxies 
they analyze. For each galaxy we then identify the low- 
est pressure Pmin for which all the pixels / rings have 
i?H2 > ^H2,min- At pressurcs above Pmin, the sample 
should be nearly complete, in the sense that no pixels 
or rings will have molecular non-detections that could 
bias a fit of Ph2 versus P. At pressures below Pmin 
the data are incomplete, and we therefore drop all pix- 
els or rings whose pressures are below Pmin- This pro- 
duces a sensiti vity-co rrected set of synthetic data. Fi- 
nally, following IBR06I we then bin the pixels by pressure, 
using bins 0.1 dex wide, and for each pressure bin in each 
galaxy we compute an a verage value of Ph2 • We do not 
bin the rings in the IL08I sample. 

We overplot the synthetic and real data in Figure O 
Fitting our synthetic data to a powerlaw function for 
Ph2 versus P gives 

i?H,= f ^ (54) 

V9.0 X 103 K cm -y ^ ^ 

if we include our models for both the lBROGi and lLOSi data 

sets. Using only one or the other gives 

V2. 8x105 K cm ' ^ ' 

As expected, we obtain a slope shallower than unity as 
a result of the systematic increase of p* with P. The 
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Fig. 9. — Molecular to atomic ratio i?H2 versu s pressure P. We 
show the data sets of iBROQ (squares) and LOS' (diamonds), and 
predictions from theory for the value of Eg and Z' fo r each data 
point (plus signs). In the top panel we show only the IBROSi data 
and the theoretical predictions corresponding to it. In the middle 
we show only the|L08|data and predictions, and in the bottom panel 
we show both data sets together. In all panels we only plot rings 
or points for which both molecular and atomic gas are detected for 
the observations, and the theoretical points are corrected for finite 
sensitivity as described in the text. We do not plot theoretical 
predictions for galaxies without measured metallicities^ We also 
show the best powerlaw fit s of BR06 data (equation 1501 dashed 
line), the LOS data (eguation lSlI dot-dashed line), and to the model 
predictions (equation l54| solid line). In the upper panel the th eory 
line uses only the B ROq data, in the middle it uses only the ILOSl 
data, and in the bottom panel we show a line fit to our predictions 
for both data sets. The dotted horizontal lines show ijjjj = 1; 
below this line our predictions should be taken as upper limits 
only 

scatter in the real and synthetic data are also compara- 
ble. As the Figure shows, our model predictions overlap 
wit h the o bserved data reasonably well, particularly for 
the IBR06I data; the best fit for the IBR06I data set is rea- 
sonably close to the observational best fit, while for ILOSl 
the theoretical predictions and the best fit to them give 
a slope similar to the observed value, but are systemat- 
ically shifted to lower i?H2- The displacement between 
our model and the'LOS' data probably exists for the same 
reason that our predictions for ^-nd R^^ in § 14.1. 2! 
are lower limits: spatial averaging leads to an underesti- 
mate of the true column densities of the atomic- molecular 
complexes in a galaxy, which in turn leads us to slightly 
underpredict the molecular content. The averaging is 
significantly worse if done over rings than over individ- 
ual beam pointings, so the offset is noticeably larger for 
the ring-averaged data set. 
It is also worth cautioning that in performing this 



fit we have included model points where we predict 
i?H2 < 1, a region of parameter space where the pre- 
dictions of our theory should be taken as upper limits. 
If we exclude these points from both our model predic- 
tions and from the observed data, all the slopes become 
significantly more shallow, but they remain consistent 
with one another. In this case our model predictions 
give a best fit of Rh^ = [(P/fcB)/(83 K cm~'^)1°- 4^ for 
the combined BR06 and L08 samples, while the IBR06I 
and [LOI data sets give [(P/fcB)/(120 K cm-^)]0-5i 
[(P/fcB)/(300 K cm-^)]0-53, respectively. These fits are 
performed without weighting by the errors or properly 
including the effects of upper limits, so they should only 
be taken as general indications, but the results do show 
that even if we limit our fit to the part of parameter space 
where we can apply our theory with high confidence, we 
obtain good agreement between model predictions and 
observations. 

We can further check the level of agreement between 
the model and the data, and see how well o ur theo ry com- 
pares to the purely empirical fit, using the IBR06I sample, 
where there is no (or less) systematic offset due to unre- 
solved clumping. In Figure [TOl we plot the same data as 
in Figure [51 but show the data and model predictions for 
the galaxies one by one. As the plots show, our model not 
only fits the general trend between i?H2 and P, for most 
galaxies we obtain a good match on a bin- by-bin basis. 
Even in those galaxies where the agreement is poor, our 
model still gives a correct lower limit for i?H2 • The level 
of agreement between our predicted curves and the data 
is as good as the purely empirical fit between i?H2 and 
pressure in IBROGI . without the need for any free param- 
eters. It is not clear why one model or the other does 
better for particular galaxies. Our model appears to pro- 
vide a substantially better fit for NGC3627, NGC5033, 
and NGC5055, while doing noticeably worse for IC 10, 
NGC0 598, and NGC4736; for the rest of the sample the 
IBR06I empirical fit and our model calculations are nearly 
equally good fits. It is not clear what galaxy properties 
favor one model or the other. 

4.1.4. The Small Magellanic Cloud 

As noted ab ove, the SM C's metallicity of 12-1- 
log(0/H) = 8.0 (|Dufouilll98"i l makes it the lowest metal- 
licity galaxy in our sample for which we have detections 
rather than upper limits on the II2 content. To analyze 
this data set , we take the maps of Hi and II2 column de- 
termined by iLerov et"aLl (|2007l ) and break the data into 
20 bins in total gas column (including He). In each bin 
we compute the mean Hi and H2 column density. It is 
worth noting that the technique used to infer H2 column 
densities in this data set have ^ 50% systematic uncer- 
tainties, so we must proceed with this caution in mind. 

To test our theoretical model, we generate a prediction 
for the H2 mass fraction as a function of total gas column 
density Scomp as in equation (j48[) , using a metallicity rel- 
ative to Solar of Z' « 0.2. Before comparing the model 
curve to the data, we face the difficult problem of incli- 
nation correction. The SMC is a triaxial structure with 
an aspect ratio of roughly 1:2:4, with 4 representin g 
the direction along the line of sight (|Crowl et al.ll200ll) . 
The large extent along the line of sight means that the 
inclination correction is very significant, but the triaxial- 
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Fig. 10. — Same as Figure[9] except that we show only the IBR03 , and we plot only two galaxies per panel, as indicated, so that we can 
compare on a galaxy- by-galaxy basis. I n each panel the black open points are the observed data, while the gray filled points are our model 
predictions. The dashed line shows the IBR03 empirical fit. 



ity of the galaxy means that no single number describes 
the inclination, as is the case for a disk. We therefore 
compare the model and the data using two different in- 
clination corrections, one assuming an inclination i — 76° 
(corresponding to a 1 : 4 aspect ratio) and one assuming 
an inclination i = 63° (corresponding to a 1 : 2 aspect 
ratio). These should bracket the true inclination cor- 
rection. Another caution we should make here is that, 
due to these large inclinations, there may be significant 
amounts of gas along a given line of sight that are not 
associated with whatever atomic-molecular complex it 
intersects, so the Hi column densities may be overesti- 
mates (see S I4.1.2|) . 

Figure [Tl] shows a comparison of the data and the 
model curve. As expected, the data for the two differ- 
ent inclination corrections bracket the model curve quite 
well. We should be cautious about reading too much into 
the agreement at low Eobs ( ^ 30 Mq pc~^), where both 
the data and our model are quite uncertain. Indeed, i?H2 
is such a sharp function of Sobs at low column that the 



factor of 2 variation in the total column introduced by 
the uncertainty in the inclination translates to a factor 
of ^ 50 difference in the value of i?H2 at a fixed Sobs- 
Thus, our model need only be correct within a factor of 
50 to lie in between the two curves for the different incli- 
nation corrections. At higher Sobs ( ^ 80 Mq pc~^), on 
the other hand, the inclination correction only produces 
a factor of ~ 2 uncertainty in i?H2 ? and the ability of our 
model to match the data here is significant. It shows that 
we are capable of predicting the total column at which 
i?H2 1 to better than a factor of 2 accuracy even in a 
galaxy with metallicity 1/10 Solar. 

4.2. Galactic Observations 
4.2.1. Data Sets 

We next compare our model to observations of the 
molecular content of clouds i n the Milky Way, as 
measured by the Copernicus (Spitz er fc Jenki"iii Il975| ) 
and the Far Ultraviolet Spectroscopic Explorer (FUSE) 
(jMoos et all 120001 : ISahnow et al.ll2000f) missions. These 
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Fig. 11. — Column densities of Hi and H2 {upper two panels) and 
H2 to Hi ratio {bottom panel) as a function of total column density 
^obs in the SMC. The solid curve shows the model prediction, 
while the dashed lines show our confidence regions, as in Figures 
IHf6l The asterisks and plus signs represent the observed column 
densities using inclination corrections of 63° and 76°, respectively. 

satellites measured absorption of ultraviolet light from 
background stars, and in some cases AGN, in the Lyman- 
Werner bands, enabling them to estimate the population 
of molecul ar hy drogen along a given line of sight. 

Savage et alj (|l977) and Bohlin ct al. ( 1978) (for 
Copernicus) and iRachford et all (|2002l 12008) (for 
FUSE) report measurements of the molecular hydro- 
gen columns in the Milky Way disk for lines of sight at 
low galactic latitude (5 ^ 10°), using stars as background 
sources. They combine these with Hi observations along 
the same lines of sight to determine mo lecular hydro- 
gen fr actions for clouds these sighthnes. iGillmon et al.l 
(|2006f ) report FUSE observations of sightlines at high 

galactic latitude (6^20°) using AGN as background 
sources. These observations probe clouds lying above 
or below the galactic plane, which are presumably illu- 
minated only on one side by st ars in the disk. As with 
the low-latitude samples, Gill mon et al.l combine these 
observations with measurements of the Hi column along 
the same sightlines to determine molecular fractions. 

These observations are different than the extragalactic 
ones in that they are pencil-beam measurements rather 
than averages over atomic-molecular complexes. We pos- 
sess little information about the geometry of the clouds 
these beams probe, and this ignorance complicates com- 
parison of the data to our model. A given observed total 
(atomic plus molecular) column density Sobs might be 
the result of a line of sight passing through the center 
of a small cloud, or might be the result of a beam that 
passes only tangcntially through a much larger cloud. 
Our models do not predict the same H2 content in these 



two cases, so in the absence of additional information 
there is no way to map a given Eobs to a unique predic- 
tion for the molecular fraction. Similarly, there is good 
evidence in the data sets from the rotational excitation 
of the H2 that at least some lines of sigh t probe multiple 
cloud s that are separated in space (e.g. iBrowning et al.l 
l2003f ). Our model predicts a lower molecular fraction 
for this case than for the case of a single cloud with the 
same Eobs- Attempts to remove these effects statistically 
are complicated by the fact that the FUSE and Coper- 
nicus lines of sight are not an unbiased sample. Each 
of the samples we use was chosen specifically to probe a 
certain range in column density or other properties, and 
all the selections are biased against high column densi- 
ties and high molecular fractions, both of which produce 
high extinctions of LW photons that make determining 
a molecular column very costly or altogether impossible. 
Below we discuss how we deal with the problems of ge- 
ometric uncertainty and selection bias for the galactic 
plane and high-latitude samples. 

A second complication for this data set is that the as- 
sumption we make in our model that clouds are subjected 
to a rela tively uniform backgro und dissociating radiation 
field fcf. iKrumholz et ani2008D may not be appropriate 
in this context. Treating the ISRF as uniform is rea- 
sonable for the giant cloud complexes with masses > 10"* 
Mq and sizes > 100 pc probed by extragalactic measure- 
ments; these sample the radiative output of many stars 
and star clusters. In contra st, the entire sample of high 
latitude clouds observed by IGillmon et al.l (|2006l ) has a 
total jnass ~ 3000 A£n spread out over a ~ 100 pc^ area 
(jGillmon fc Shul]||2006f) . making them tiny m compari- 
son to giant complexes. The galactic clouds probed by 
the Copernicus and FUSE observations are of unknown 
size, but are likely to be similarly small. This is because 
any line of sight passing through a true giant cloud com- 
plex would be completely opaque in the Lyman- Werner 
bands, thus preventing Copernicus or FUSE from mak- 
ing any measurement of the H2 column. The dissoci- 
ating flux in the vicinity of such small clouds may be 
dominated by a single nearby star or star cluster. 

Due to the effects of uncertain geometry and non- 
uniform radiation fields, we expect this data set to show 
significantly more scatter than the extragalactic one. 

4.2.2. The Galactic Plane Data 

Given the limitations of pencil-beam measurements, 
we must compare our model to the data in a way that ac- 
counts for our geometric uncertainty. The most straight- 
forward comparison we can make is by considering a lim- 
iting case. Consider an observation that measures a to- 
tal column density Eobs- According to our model, the 
fraction of gas in molecular form will depend on the ge- 
ometry of the complex; if the line of sight probes the 
outer parts of a very large cloud, always staying within 
the atomic region, then in principle there could be no 
molecules present no matter how large Sobs might be. 
The reverse is not true, however. Even if the line of 
sight passes exactly through the center of an atomic- 
molecular complex, it must still pass through atomic 
shielding layer, and there will therefore be some Hi along 
that line of sight. In our model, therefore, there is a min- 
imum amount of Hi that we predict must be present for 
a given column of H2; this minimum corresponds to the 
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Fig. 12. — Column density ol H2 {upper panel) and Hi {lower 
panel) versus total column density for the Copernicus {diamonds) 
and FUSE {squares) data, and for the limits {dashed line) and 
typical values {solid line) we compute from our model. The dashed 
line is an upper limit in the upper panel, and a lower limit in the 
lower one. For the data, we only plot points with measured Hi 
columns. Arrows indicate lines of sight for which the Hi column is 
measured but only an upper limit is available for the H2 column. 
In the upper panel the dotted lines show Sjjj = S and Sjja = 
S/3 (the edge of our region of confident prediction), while in the 
lower panel they show Ejji = S and Ejji = 2E/3. Note that Afjjj 
represents the column of H nuclei in molecular form; the number 
of H2 molecules is half this. 



case of a line of sight that passes directly through the 
center of a complex. This implies that we can make a 
prediction as follows: we consider a complex of mean 
column density Ecomp, and use our model to predict the 
atomic, molecular, and total column densities, Shi, con, 
5^H2,cen: and Scon: that One would see along a line of sight 
passing through the center of the complex. The curves 
Sni.cen versus Seen and SH2,con versus Seen that we gen- 
erate through this procedure should then be lower and 
upper limits, respectively, on the observed distributions 
of Shi versus Sobs and Sh2 versus Sobs- 

Instead of computing a limit, we can also make an esti- 
mate for the "typical" atomic and molecular fractions we 
should see at a given Sobs- If we knew the true distribu- 
tion of column densities for atomic-molecular complexes 
in the Milky Way, we could do this by integrating all 
lines of sight through that distribution and computing 
the mean atomic and molecular columns at a given total 
column. However, we do not know the true distribution, 
and even if we did this procedure would not address ef- 
fects of observational bias in selecting sightlines. Instead, 
we make a much rougher calculation. Using our model we 
can compute the radius XHaCSeomp) at which the cloud 
whose true mean column density is Seomp transitions 
from atomic to molecular for solar metallicity {Z = 1). 
The atomic and molecular columns along a sightline that 
strikes the cloud of radius i? at a distance /3R from its 
center are 
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Fig. 13. — versus column density for the FUSE high-latitude 
sample ( .sqMares; G illmon et al..,200S) and for our model {line). Ar- 
rows indicate lines of sight for which only an upper limit is available 
for the H2 column. Note that N^hs represents the total column of 
H nuclei in either atomic or molecular form. 
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and the II2 fraction averaged over all pencil beams pass- 
ing through the cloud is, 



/H2 



jbcams (Seomp 



-'HI 



-df3. 



(58) 



To estimate a "typical" atomic or molecular content 
for a given observed total gas column Sobs, we can 
simply take SnLobs « [1 - /h2, beams (Sobs)] Sobs and 
SH2,obs ~ /H2,boams(Sobs)Sobs- This is equivalent to 
saying that our sightlines do indeed probe random im- 
pact parameters on cloud complexes, and that when 
we observe a sightline of column density Sobs, most of 
the time we are observing a complex whose mean col- 
umn density is also about Sobs- This assumption could 
fail if low column density clouds were rare compared to 
high column density ones, so that most of the sightlines 
that produce a given Sobs were really tangential paths 
through high column clouds rather than beams passing 
close to the the centers of low column clouds. However, 
there is no evidence that high-column density atomic- 
molecular complexes outnumber low-column ones, and 
in general in the ISM more mass tends to be in diffuse 
than dense structures, so we proceed with our assump- 
tion that Seomp ~ Sobs- 

We plot both our limits on the Hi and H2 columns and 
our estimates for their typical values against the Coper- 
nicus and FUSE data sets in Figure [121 As the plot 
shows, both our limits and our typical columns match 
the observed data reasonably well. As with the extra- 
galactic case, we recover the overall trend that the gas is 
mostly Hi until a total column of a few Mq pc~^, and 
then mostly molecular thereafter. 

4.2.3. The High-Latitude Data 

The comparison to the high-latitude data set is some- 
what less straightforward, because the clouds at high 
latitude are o n ly illu minated from one side. Indeed, 
iGillmon et al.l ()2006f ) find that these clouds become 
molecular at column densities a factor of ~ 2 lower than 
do clouds in the galactic plane, probably as a result of 
this one-sided illumination. We have not considered the 
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case of clouds of finite size illuminated over only half 
their surfaces. However, we can obtain a reasonable ap- 
proximation to this configuration using the case of semi- 
infinite clouds subjected to an isotropic, uniform disso- 
ciation radiation field, which we considered in § [2l A 
semi-infinite cloud blocks dissociating radiation over 2tt 
ST, and thus the depth of the atomic shielding layer at its 
surface is determined only by photons that arrive from 
the "front side" . This is therefore a close analogy to the 
case of a high-latitude cloud, although it differs in the 
cloud geometry. 

In Paper I we show that using a semi-infinite slab in 
place of a finite cloud tends to produce errors in esti- 
mating the exact location of a transition from atomic 
to molecular in a finite cloud, but that the semi-infinite 
calculation does give a good estimate of the column at 
which clouds transition from mostly atomic to mostly 
molecular. We therefore expect to obtain roughly the 
right transition column and thus the right column of Hi 
for clouds that are mostly molecular, but not particularly 
accurate predictions for the exact molecular column. For 
a semi-infinite cloud of solar metallicity, as we show in 
§ [2] we predict a molecular to atomic ratio of 

n _ / 0, S < 4.5 Mq pc-2 

~ \ (S-4.5Af0 pc-2)/l], S>4.5Mopc-2 ' 

(59) 

where S is the total column density. In other words, 
for a semi-infinite cloud of solar metallicity, the first 4.5 
Mq pc~^ are atomic, and the rest are molecular. We 
show this prediction overplotted with the FUSE high- 
latititude data in Figure[T51 As the plot shows, we indeed 
do not get very accurate predictions for the exact ratio of 
Hi to H2, but our calculation agrees quite well with the 
general value of S for which the transition from molecu- 
lar to atomic occurs. In particular, our one-sided model 
recovers the observational result that a cloud illuminated 
from one side shows a lower atomic-to-molecular transi- 
tion column than an isotropically-illuminated cloud. 

5. SUMMARY AND CONCLUSIONS 

We present a first-principles calculation of the molec- 
ular gas content of galactic disks in terms of the observ- 
able properties of those galaxies. Our calculations build 
on the simple model for photodissociation fronts in finite 
clouds presented in Paper I, in which we show that the 
amount of atomic material required to shield a molecu- 
lar cloud against dissociation by the interstellar radiation 
field (ISRF) can be characterized by two parameters: Xj 
a radiation field strength normalized by the gas density 
and the properties of dust grains, and tr, a measure of 
the dust optical depth of a cloud. We show that, due to 
the way the density in the cold phase of the atomic ISM 
varies with ISRF, the normalized radiation field strength 
takes a characteristic value % « I in all galaxies where 
a two-phase atomic ISM is present, with only a weak 
dependence on metallicity. 

The existence of a characteristic normalized radiation 
field strength, and its weak dependence on metallicity, 
has a number of important consequences. First, it en- 
ables us to give a simple analytic approximation (equa- 
tions [36] or [37]) for the fraction of mass in an atomic- 
molecular complex that will be in the atomic or molec- 
ular phases solely in terms of the column density of the 



complex and the metallicity of the gas. This makes it 
easy to test our calculations against observations. Sec- 
ond, we show that as a consequence of x assuming a 
nearly fixed value, the atomic envelopes of molecular 
clouds have a characteristic visual extinction Ay ~ 0.2 
at solar metallicity. Similarly, the transition from atomic 
to molecular gas occurs at a characteristic shielding col- 
umn of Shi ~ 10 M© pc~^. These quantities both vary 
sub-linearly metallicity, with Ay declining and Shi in- 
creasing as metallicity does. We calculate these values 
and their metallicity-dependence solely in terms of the 
microphysical constants that describe the properties of 
H2 formation and dissociation, the cooling curves of Cii 
and Oi, and the properties of interstellar dust grains. 
Our model does not depend on unobservable parameters 
such as the gas volume density or the ISRF strength in 
a galaxy. 

Our model compares favorably with observations of 
the atomic and molecular content of clouds both in the 
Milky Way and in nearby galaxies. We are able to re- 
produce both the characteristic column density at which 
clouds transition from being primarily atomic to primar- 
ily molecular, and the way that this characteristic column 
depends on metallicity and on whether clouds are illumi- 
nated on one side or on both sides. We are also able 
to reproduce the observed correlation between molecular 
content and interstellar pressure. 

The development of a predictive model for the molecu- 
lar content of galaxies that does not rely on unknown and 
generally unobservable quantities such as the radiation 
fields or gas volume densities in those galaxies opens up 
new possibilities to advance our understanding of galactic 
evolution. In low-density dwarfs or high-redshift galaxies 
containing few metals, the formation of molecular clouds 
may be the rate limiting step in star formation. On the 
other hand, previous work has shown that, once molec- 
ular gas forms, it converts itself into stars at a rate of 
a few percent of the mass per free-fall time ind ependent 
of its density or environment (iKrumholz fc McK cc 200^ 
IKrumholz fc"TMl [2007l:[McKee fc Ostrikeill2007l) . Thus 
a theory of molecule formation creates the possibility of 
developing a theory of the star formation rate capable of 
making predictions that can be applied not only in the 
relatively molecule-rich nearby galaxies for which em- 

f'rical star formation laws have been determined (e.g. 
ennicutt 1998), but also in the more distant and lower 
metallicity uni verse where thes e laws are known to break 
down (e.g. Wolfe fc Chenll2006f ) . We plan to explore such 
a theory in future work. 
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